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Abstract Our knowledge of the high-energy universe is 
undergoing a period of rapid change as new astronom- 
ical detectors of high-energy radiation start to operate 
at their design sensitivities. Now is a boomtime for high- 
energy astrophysics, with new discoveries from Swift and 
HESS, results from MAGIC and VERITAS starting to 
be reported, the upcoming launches of the 7-ray space 
telescopes GLAST and AGILE, and anticipated data re- 
leases from IceCube and Auger. 

A formalism for calculating statistical properties of 
cosmological 7-ray sources is presented. Application is 
made to model calculations of the statistical distribu- 
tions of 7-ray and neutrino emission from (i) beamed 
sources, specifically, long-duration GRBs, blazars, and 
extagalactic microquasars, and (ii) unbeamed sources, 
including normal galaxies, starburst galaxies and clus- 
ters. Expressions for the integrated intensities of faint 
beamed and unbeamed high-energy radiation sources are 
also derived. A toy model for the background intensity 
of radiation from dark-matter annihilation taking place 
in the early universe is constructed. Estimates for the 7- 
ray fluxes of local group galaxies, starburst, and infrared 
luminous galaxies are briefly reviewed. 

Because the brightest extragalactic 7-ray sources are 
flaring sources, and these are the best targets for sources 
of PeV - EeV neutrinos and ultra-high energy cosmic 
rays, rapidly slewing all-sky telescopes like MAGIC and 
an all-sky 7-ray observatory beyond Milagro will be cru- 
cial for optimal science return in the multi-messenger 
age. 
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1 Introduction 

The next decade is likely to be remembered as the pio- 
neering epoch when the first high-energy (PeV - EeV) 
v sources were detected with IceCube pQ and its km- 
scale Northerm hemisphere counterpart, and when the 
problem of cosmic-ray origin was finally solved through 
identification of the sources of cosmic rays at all ener- 
gies, from GeV - TeV nucleonic cosmic rays accelerated 
by supernova remnant shocks of various types, to extra- 
galactic super-GZK 7-ray and v sources. 

The cosmology of 7-ray sources in the s» 10 MeV - 
10 GeV range is treated here. The lower bound of this 
energy range ensures that the 7 rays originate from non- 
thermal processes, and the upper bound is defined by 
the energies of photons that originate from sources at 
redshifts z> 1 without significant 77 — » e + e~ attenua- 
tion in reactions with photons of the extragalactic back- 
ground light (EBL). The formalism also applies to other 
nonthermal radiations, to ultra-relativistic particles, in- 
cluding PeV - EeV v and ultra- high energy neutrals, 
and to multi-GeV - EeV photons by taking into account 
attenuation and reprocessing of the 7-rays on the EBL. 

The problems treated here are the 

1. Event rate of bursting sources; 

2. Size distribution of bursting sources; and 

3. Apparently diffuse intenstity from unresolved sources. 

I outline applications of these results to beamed souces, 
including GRBs, blazars and extragalactic microquasars, 
and unbeamed sources, including star- forming galaxies 
and merging clusters of galaxies. 

This paper, prepared for the conference proceedings 
of the Multi-Messenger Approach to High Energy Gamma- 
Ray Sources, held 4-7 July 2006 in Barcelona, Spain, 
addresses in a more formal manner the points I was to 
cover, including blazars which I could not neglect (see 
Ref. [2] for a review of blazar emissions). The formalism 
applies to analysis of 7-ray and v data from GLAST, Ice- 
Cube, and other high-energy astroparticle observatories. 



2 



Charles D. Dermer 



2 Event Rate of Bursting Sources 

The Robertson- Walker metric for a homogeneous, isotro- 
pic universe can be written as 



ds 2 



c 2 dt 2 



dr 2 



kr 2 



'-dQ) 



(1) 



where r is a comoving coordinate and R(t) is the ex- 
pansion scale factor. The most convenient choice is to 
have r take the value of physical distance at the present 
epoch so that R(t) = R = 1, and denote i?» = R(t*) at 
emission time £* < t (stars denote the emission epoch). 
Material structures reside for the most part on constant 
values of the comoving coordinates, whereas light and 
ultra-relativistic particles cannot be confined to such co- 
ordinates. From the definition of redshift z = (A— A*)/A*, 
we have 1 + z = e*/e = At/ At* = where e refers 

to the energy of the photon or ultrarelativistic particle. 
The curvature of space is determined by the curvature 
constant k, with k = for flat space. 

The proper volume element of a slice of the universe 
at time t» is, from eq. ([T]) for a flat universe, 

dV* = Rlr 2 dr = dr^dQ* = cdt*dA* . (2) 

Comparing with the definition dA — (Rr) 2 dQ and noting 
that dQ* = dQ in the absences of cosmic shear, we have 



dA, 
~d~A 



1 



(l + z) 2 ■ 

The directional event rate, or event rate per sr, is 
dN 1 f . dt 



(3) 



dz 



. dt* (R*r) 2 h*(z) 



dz 



(1 + *) 



(4) 



where the burst emissivity n*(z*) gives the rate density 
of events at redshift z. An expression for (R*r) 2 can be 
derived by recalling the relationship between energy flux 
<I>E and luminosity distance d^, namely 



d£ 



® E 4ird 2 ^d 2 dU 



dAdt 

so that with eq. ((3j) 



(R*rf 



d 2 L (z) 



For a flat vlCDM universe, 
dz 



(I + z) 2 



(6) 



|— | = h (i + z)y/n m (i + z)3 + n A 



(7) 



[5P|. where H a = 72 km s" 1 Mpc" 1 , Q m = 0.27 and 
Qa = 0.73 are the ratios of the energy densities of total 
mass, including both normal matter and dark matter, 
and dark energy, respectively, compared to the critical 
density for the flat /1CDM cosmology of our universe. 
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Fig. 1 Different structure formation histories (SFHs) consid- 
ered in this paper. As labeled, CCSFH: constant comoving 
SFH; LSFR: lower star formation rate; USFR: upper SFR; 
SFR HB: SFR history from Ref. [7]; SFH IR: SFH of IR lu- 
minous galaxies 11 ; GRB SFHs: range of SFHs of GRBs 
used to fit Swift and pre-Swift GRB distributions [6]; SFH 
BL: SFH of BL Lac objects 0. These rates are poorly known 
at z > 1. 
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The directional event rate, eq. dU), becomes 

dN_ _ f°° ^dt^ d 2 L (z)h co {z) 
dQ C J Q l dz l (l + z) 2 ' 

after using the relation n* — (1 + z) 3 n co (z) to write the 
directional event rate in comoving rather than proper 
quantities. If separability between the emission prop- 
erties and the rate density of sources can be assumed 
(a crucial assumption), then n co ^(z) = hiUi(z), where 
Si(z) is the structure formation history (SFH) of sources 
of type i, defined so that Si(z — > 0) = 1, and hi is the 
local (z <C 1) rate density of bursting sources of type i 
(see Fig. 1). 

The comoving density can be formally expanded as 



hco(z) 



dQ' 



daN(a;z) / . . . h co {Q' , a, . . . ; z)(9) 



0[6]. The direction Q 1 = (p' = cos #',</>') specifies the 
orientation of the system with respect to the direction 
to the observer, and N(a; z) is a normalized distribution 
function for parameter a. For example, a could represent 
the bulk Lorentz factor r, the total energy radiated, the 
comoving- frame power, or the spectral index of the radia- 
tion. For sources oriented at random, h co (Q' , cv, . . . ; z) — 
h co (a, . . . ; z)/4tt. For two-sided jet sources, —1 < fjf < 1 
and h co {f}' , a,...;z) = 2h co (a, . . . ; z)/4tt. 

If one were to consider persistent rather than burst- 
ing sources, an analogous derivation gives the directional 
number count of sources of the type i as 
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3 Size Distribution 

After substituting eq. © into eq. © and placing limits 
on the integrals in accordance with detector specifica- 
tions, model distributions of source properties can be 
derived, in particular, the size distribution. 



3.1 Beamed Sources 

A distinction between two types of models for relativis- 
tically beamed sources needs to be made. A blast-wave 
model is usually considered for GRB sources. Here the 
half-angular extent 9j of the collimated spherical blast- 
wave jet is assumed to be much greater than the Dopplcr 
beaming angle 9b ~ so that 9j ^> 9b, where r is 
the bulk Lorentz factor of the outflow. In this case, the 
observer is limited to detection of a GRB if the direction 
from the source to the observer intercepts the solid an- 
gle of the blast wave. By contrast, in a blob model, as is 
usually considered in blazar studies, 6j <C 9b, and the 
beaming properties of the jet are determined primarily 
by the Dopper factor 



[r(i-/3 M )]- 



(ii) 



where (3T = \JT 2 — 1, and arccos /i is the observer's an- 
gle measured with respect to the jet axis. Interpreting 
emissions within a blast-wave and blob framework pro- 
vide the simplest models that can be used to system- 
atically analyze the statistics of GRBs, blazars, radio 
galaxies, and microquasars within the context of a phys- 
ical (rather than a phenomenological) model. 

Gamma Ray Bursts If a GRB releases an amount of 7- 
ray energy £* 7 that is deposited in a waveband to which 
a GRB detector is sensitive, then an event is recorded 
when the source flux 



47rd|(z)(l - fj,j)At*\ b 



(12) 



where f e is the vF v threshold sensitivity of the GRB de- 
tector. Beaming of the emission into a jet with opening 
half-angle 9j has the effect, for constant £„ 7 , to enhance 
the received flux by a factor (1— though the chance 
of the jet being in the line of sight to the observer (com- 
pared to an isotropically emitting source) is reduced by 
the factor 1 — /%•. The term At, is a bolometric correction 
factor made in lieu of a full spectral treatment. 

The GRB size distribution for the blast-wave geom- 
etry is given by 



dN GRB {> U 
dfl 



2cn 



GRB 



dz 



dU , d 2 L (z)S GRB (z) 
1 dz 1 (1 + z) 2 



x / dfj,j gQij) (1 - Hj), 

J max(0,/ij ) 



(13) 



where g(nj) is the normalized distribution function of 
GRB jet opening angles, and 
f 



l l J 



1 



4nd 2 L (z)At*\ b f e 



(14) 



Le and Dermer (2006) [5] have used this approach 
to analyze the redshift and opening-angle (9j) distribu- 
tions of GRB detectors, including missions before Swift 
compared with distributions measured with Swift. They 
find that the comoving rate densities of GRBs must un- 
dergo positive evolution to at least z ^ 5 - 7 to account 
for the difference in distributions of pre-Swift and Swift- 
detected GRBs with redshift information. By contrast, 
the star-formation history of the universe as inferred 
from blue and UV luminosity density, peaks at z ~ 2 - 3 
and seems to decline at larger redshift [7] . Le and Dermer 
(2006) find this SFH to be incompatible with the statis- 
tics of GRBs with measured redshifts. Thus the SFH of 
GRBs is apparently very different than the integrated 
high-mass star formation history of the universe. 

This approach can be suitably adapted to the short, 
hard class of GRBs to infer the rate density of this class 
of GRBs. A large data set, accumulated after a long 
Swift lifetime, can in principle distinguish between mod- 
els involving compact-object coalescence and accretion- 
induced collapse of neutron stars. 

Blazars A considerable simplification to the emission prop- 
erties of blazars results by approximating the vF u fluxes 
detected from a distant source by the expression 



> fe 



(15) 



rproc e D^z 

h ~ d 2 L (z) 

[S11S], where l' e is the comoving directional power and 
f e represents the characteristic flare size, in this case, 
in units of energy flux. The beaming factor indices for 
individual radiating blobs are 

f (p + 5)/2, synchrotron/SSC 



.P+'A 



EC 



The synchrotron/SSC factor applies to blazars where the 
7 rays are from the synchrotron self-Compton processes, 
specifically X-ray-selected blazars and TeV blazars. The 
external Compton (EC) beaming factor applies to blazars 
where the 7-rays are ambient photons, external to the jet, 
that intercept the jet and Compton-scattered by the jet 
electrons. Examples of ambient radiation fields are the 
accretion disk photons, and accretion-disk photons that 
are scattered by surrounding dust and gas. 

From eqs. © and ^j, the blazar flare size distribu- 
tion is given by the expression 



dN b , 
df2 



■(>/«) = 2cn u I , h ^d 2 L {z)E bl {z) 



dz' (1 + z) 2 

/oc poo 
dTN(r-,z) J d£' e N(£' e ;z) [1 - max(-l, /*)] ,(16) 
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where 

*=^-kW° ] - (17) 

Specification of the z-evolution of the normalized dis- 
tribution functions N(r; z) and N(£' e ; z) due to number 
evolution or luminosity evolution, respectively, connects 
this formulation back to the cosmology of physical pro- 
cesses and the growth of structure taking place in the 
early universe. 

Refs. [5] uses this approach to analyze the redshift 
and size distribution of EGRET 7-ray blazars (see also 
[5]), divided into flat spectrum radio quasars (FSRQs) 
and BL Lac objects (BLs). Evolutionary behaviors are 
found that characterize the measured redshift and size 
distributions of FSRQs and BLs. The behavior of the BLs 
is in accord with the conjecture that BLs are late stages 
of the formation and evolutionary history of FSRQs, and 
before that, IR luminous galaxies [TUlfTTj . See Ref. [5] for 
predictions of the number of blazars that GLAST will 
detect. 



Microquasars 7-ray emission from microquasars could 
be visible from nearby galaxies if the bulk Lorentz fac- 
tors in microquasar jets were large enough that the re- 
ceived flux from a microquasar in another galaxy was 
brighter than threshold. The size distribution of micro- 
quasar flares can be written by taking the limit z <C 1 of 
the blazar expression, eq. (|f 6[) . to give 



and heavily extincted starburst and infrared luminous 
galaxies (whose ambient radiation would attenuate the 
7 rays), we can count the number of source detections 
above a threshold flux / e , following eq. (fTT)]) . to give: 



dNj 



(>/e) 



00 dKz^jz) 
dz 1 (l + z) 



dL* N(L*;z) 



(20) 



The luminosity function of the unbeamed source popu- 
lation is denoted by N(L*;z), where = J Q de*Z/(e*) 
is the total luminosity of the source. Writing the spec- 
tral luminosity L*(e*) = L* e~ 1+a " gives the threshold 
condition 



And 2 ~ h 



(21) 



for detection of these sources. For a power-law spectrum 
with low- and high-energy cutoffs, this expression can be 
used to impose the lower limit L„ nin in eq. (|20[) . which 
also assumes an average over large volumes. For normal 
galaxies, volumes of radii of several Mpc may be large 
enough for this averaging. For clusters of galaxies, an 
averaging size scale of many tens of Mpc is needed, as 
calculations at scales less than z ~ 0.02 are subject to 
strong fluctuations due to the low density of clusters of 
galaxies in this volume. 



da 



■(> /«) 



2c 3 n M 
n 



dz z 2 S m (z) 



dr N(r; z) [1 - max(-l, /})] 



(18) 



which assumes an averaging over the small scale mass 
distributions of nearby galaxies. The integration in l' e 
is removed in this expression, compared with eq. (fTB|) . 
by assuming an Eddington limitation on the accretion 
flow. Approximating the emission spectrum by a single 
power law with vF v index a v in the comoving energy 
range e' < e' < e[ , the directional luminosity is therefore 
limited by 



pi < 



2 x 10 d8 m c 



ergs s sr 



(19) 



noting that the emission is beamed into ps r~ 2 of the 
full sky, and that the radiated power is boosted by r 2 
due to bulk motion of the plasma. Here iuq is the Chan- 
drasekhar mass (in units of 1.4 M Q ) of the compact ob- 
ject in the microquasar. 



3.2 Unbeamed Sources 

For 7-ray emission from unbeamed sources, like the Milky 
Way galaxy, normal galaxies, and all but the most dusty 



4 Intensity of Unresolved Sources 

The differential spectral flux 

, ,, , dN n*(e*; z)de*dt*dV* 
d0(£) = dAMde = dAMde ' (22) 

Using the relations dV* = dr^dA^ — cdt*dA/(l + z) 2 
from eq. ©, e* = e(l + z) = e z , and dt = dt*(l + z), we 
have 



(e) = c/ dz \ d ^M^ z ) 



io dz' {1 + z) 2 

Because the U i/F v " intensity el e = m e c 2 e 2 (j)(e)/ATT, 
< , _ _c_ f°° dz 1^*1 m e c 2 elh co (e^z) 



47r ,/ dz l + z 



(23) 



(24) 



GRBs The diffuse intensity of GRBs is, from eq. (J21J) 
and assuming a two-sided GRB jet source, 



m e c 3 n G RB f°° dz i dt * i s grb(z) 
4tt J dz l + z 



elf RB {< f e ) = / dz 



d Vj 5(Mi)(l - Mi) 



(25) 



For a flat vF v spectrum that covers the waveband of the 
GRB detector, m e c 2 e^A(e*; fij) = £* 7 /[Ab(l — fij)], and 
flj is given by eq. |LT 
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Using the parameters derived from analyses of statis- 
tical distributions of GRB data [6], one can then calcu- 
late the integrated 7-ray background from GRBs which, 
as we shall see, is a negligible fraction of the diffuse 
isotropic 7-ray background. Suitable scalings are adopted 
in model calculations of v -emissions from GRBs to cal- 
culate the diffuse ~ 100 TeV - EeV v intensity from 
GRBs [12] • 

Blazars The total intensity from two-sided blazar jet 
sources, which will include emission from aligned and 
misaligned blazars and radio galaxies, is given in the blob 
framework by 



2tt 



elT = r- / dz \-=- 



, dt* , 1 



dz 1 + z 



x f dn'elq b i(£*, ft'; z) 



(26) 



where qu(e*, ft 1 ', z) is the directional spectral flux of a 
blazar jet, given in the blob framework by 

e*<7bz(e* , fl'\ z) = £' e (z)nu(z)6^e"" 

[5]. Here ribi(z) is the comoving density of blazar sources. 
The intensity of unresolved blazars and radio galaxies is 
then 



el b j(<fe)- 



dz 



dt* n b i(z) 



J 1 dz 1 (1 + z) 1 



drN(r-, z) 



x / <V(^;z){[l- / 3min(l, / l)] 1 -"-(l+/3) 1 -n , (27) 



with jl given by eq. (I17p . 

Microquasars The intensity from microquasars is given 
essentially by eq. ([26]) , though with a very different local 
rate density n M<3 , SFH S^q, and distribution in r and £' e . 

Unbeamed Sources If N(r*; z) is the redshift-dependent 
luminosity function of unbeamed 7-ray sources, such as 
normal galaxies and starburst and IR luminous galaxies, 
then the diffuse intensity from these sources over cosmic 
time is 



^ I dt* J Ul e C £*^2So(^*> ^0 

dz 1 + z 



di* N(L*;z) , 



(28) 



where now L™ m (z) again depends on detector charac- 
teristics according to the prescription of eq. I|2ip. 

It is interesting to note that the factor \dt*/dz\ as- 
sociated with the passage of time in an expanding uni- 
verse saves us from Olbers' paradox. In this formulation, 
the logarithmically divergent integrated intensity emit- 
ted by radiant sources distributed uniformly throughout 
the universe is blocked by the redshifting of radiation 
and the finite age of the universe. 



GZK Neutrino Intensity The intensity of v formed as 
photopion secondaries in the interaction of UHECRs with 
the EBL is given, starting with eq. (|24|) . in the form 



el 



GZK 



„3,2 



dz 



,dt*. flGZK,*(e*, ft*i z) 



(29) 



o dz 1 (1 + z) 2 

The production spectrum of secondary v is given by 

//•oo 
dfl J de* n* h (e*, i?; z) 



dQl 



d ( M 

d 7 ; (i-co S ^)n;( 7 ;,n;-z)^f^f 



(30) 



The sum is over various channels leading to production 
of neutrinos, and n* h (e«, J7; z) and n*(j*, 0*;z) are the 
evolving EBL and UHECR proton spectra, respectively 
(generalization to ions is straightforward). Ref. [13] uses 
this formalism to calculate the GZK v intensity under 
the assumption that the sources of UHECRs are GRBs. 

The GZK 7-ray intensity can be calculated according 
to this formalism by convolving the redshift-dependent 
differential intensity with a source function that repre- 
sents the emergent 7-ray spectrum after reprocessing on 
the background radiation field. This will produce a com- 
plete model of UHECRs, which consists of a fit to the 
UHECR spectrum, a prediction for the GZK v flux, and 
the predicted diffuse 7-ray spectrum — which must be less 
than the diffuse EBL at 7-ray energies [14] . 

Dark Matter Annihilation Astrophysical searches for sig- 
natures of dark matter annihilation target regions of en- 
hanced (dark) matter density, such as cuspy cores of 
quasi-spherical galaxies, for example, dwarf ellipticals. 
Because of its proximity, even the center of the Milky 
Way is considered to be a hopeful site of dark matter 
annihilation, in spite of perturbing warps and bars in its 
normal matter distribution. One region where unavoid- 
ably high densities of dark matter had to persist was in 
the early universe. 

Up to now, we have considered classes of sources 
whose density scales as a redshift-dependent structure 
formation rate Si{z) for sources of type i. To first ap- 
proximation, source density is proportional to the total 
normal matter content, so that the factor (1 + z) 3 is re- 
moved and the SFH is described in terms of the comov- 
ing rate density. Dark matter annihilation scales as the 
square of the density of matter, so to first order we can 
write the diffuse background intensity from dark matter 
annihilation as 



el 



DM 



S 2 

m e c e 

47T 



dZl dz l (1+Z)2 ■ [6L) 

The spectral production rate of dark matter secondaries 
is written as 



Jx"dm( 1 + z) 6 a DM S(e* - e x ) 



(32) 
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where e x is the energy of secondary 7 rays or v produced 
in dark matter annihilation, and j x is the multiplicity of 
the secondaries. The maximum redshift z max represents 
the redshift where dark matter was created or fell out of 
equilibrium. 
Hence 



elP M oc 



dz 



(1 



\5/2 



-1/2 



(33) 



for e }t e x / z max . A diffuse v background from annihila- 
tion of dark matter particles with masses of ~ 10 GeV - 
TeV could peak near 10 - 100 MeV, for z max ~ 10 4 . A 
component of the diffuse extragalactic 7-ray background 
in the 10 MeV - GeV range would be formed under the 
same circumstances as the annihilation 7 rays cascaded 
to photon energies where the universe becomes transpar- 
ent to 77 pair production. 

If z max corresponds, however, to a redshift where the 
temperature of the CMB corresponds to the dark matter 
particle energy, then this indirect signature of dark mat- 
ter annihilation is probably not detectable. In any case, 
a residual photonic or v signature from dark matter an- 
nihilation in the early universe will be associated with 
any assumed dark matter annihilation cross section, and 
this emission signature cannot exceed measured values 
or upper limits. 



5 Discussion 

By inferring source densities and event rates from as- 
tronomical observations (Fig. 1), fits can be made to 
statistical (e.g., redshift and size) distributions of high- 
energy sources detected with GLAST and other high- 
energy telescopes. A model that fits the distributions en- 
tails an imperative, for each source class, to show that the 
superpositions of radiations formed by the faint model 
sources below the detection limit do not overproduce 
the measured diffuse 7-ray background and upper limits 
to the v background radiations. We illustrate the tech- 
nique in what follows, first considering detection of quasi- 
isotropic cosmic-ray induced emissions from star-forming 
galaxies. 



5.1 7-rays and v from Unbeamed Sources 

Cosmic-ray induced emissions from extragalactic sources 
are weak pl5|fl6] , which must be the case in order to agree 
with the lack of high significance detections of 7 rays 
from star-forming galaxies without jets. Other than the 
LMC 17 , which was observed with EGRET at a flux 
level consistent with that expected if cosmic rays were 
produced at a rate proportional to the star formation 
rate of the Milky Way, no unbeamed extragalactic high- 
energy radiation source has yet been detected with high 
confidence. 



Normal, Starburst, and Infrared Luminous Galaxies. 
If cosmic-ray induced emissions are primarily respon- 
sible for the high-energy quasi-omni directional emis- 
sions from extragalactic sources, then the level of emis- 
sion from the Milky Way can be appropriately scaled 
to estimate the expected flux levels of galaxies of dif- 
ferent types. The 7-ray photon production rate from 
the Milky Way inferred from COS-B observations [TH] 
is N 1 w (1.3 - 2.5) x 10 42 ph(>100 MeV) s" 1 implying 
a > 100 MeV 7-ray luminosity from the Milky Way of 
10 39 L 39 ergs s -1 , with L 3g = (0.16 - 0.32). Analysis us- 
ing GALPROP Galactic cosmic-ray propagation model 
[20] indicates that, taking into account the GeV excess in 
the diffuse galactic emission observed with EGRET and 
using a larger model Milky Way halo, the > 100 MeV 
7-ray luminosity of the Milky Way is 10 39 ergs s -1 , with 
L39 = (0.71 — 0.92). Some 90% of this emission is due to 
secondary nuclear production when cosmic rays collide 
with gas and dust in the Galaxy. 

Approximating the integrated > 100 MeV photon 
spectrum as a power law with a mean photon spectral 
index = 2.4 implies that the uL v spectrum of the Milky 
Way is 

eL M w{e) = 3.3 x 10 39 Z 39 e~°- 4 ergs s" 1 (34) 

for e = hv/m e c 2 > 200 (i.e., > 100 MeV), or N(> 100 
MeV)= 10 43 ph s _1 . 

By scaling nearby galaxies according to their super- 
nova rates, a simple estimate for the 7-ray and v emis- 
sions can be made. For example, the supernova rate of 
Andromeda (M31), at a distance of « 800 kpc, is « 1 
per century, compared to the rate of 2.5 every century 
in the Milky Way [15J . Thus the expected 7-ray photon 
flux from M31 should be at the level 

0A/3i(> 100 MeV) w 



1 10 43 ph s 



2.5 4tt(800 kpc) 2 



0.9 x 10" 8 L 39 cm" 2 s" 1 



(35) 



for a u — —0.4, which would be significantly detected 
with GLAST. Pavlidou and Fields (2001) Q2] perform 
a more detailed treatment of local group galaxies and 
predict that the > 100 MeV integral photon flux from 
the SMC, M31, and M33 are at the levels of 1.7 x 10" 8 , 
1.0 x 10~ 8 , and 0.11 x 10~ 8 ph (> 100 MeV) cm" 2 s" 1 , 
respectively. GLAST should therefore detect at least the 
SMC and M31, though its 1 year sensitivity of « 0.4 x 
10~ 8 ph(> 100 MeV) cm" 2 s" 1 means that few other 
local group galaxies are likely to be detected. The dif- 
ference between the simple-minded treatment presented 
here and their more detailed treatment is a consideration 
of the total target mass density of the different galaxies, 
and diffusion and escape of cosmic rays from the galaxy, 
which is especially important for the Magellanic Clouds. 

By extrapolating the integrated diffuse galactic con- 
tinuum emission from the Milky Way to TeV energies, 
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the integral number flux of 7 rays from a Milky- Way 
type galaxy at the distance d is 

4>{> e) S 2.4 x 10" 5 r?L 39 e" 1 - 4 /d(Mpc) 2 

M (d/1 Mpc X )^ 7 (300 GeV)]^ pM> ^ Cm " S " (36) 
where the r] factor accounts for the different supernova 
rates and target densities for the galaxy under consid- 
eration, as well as the reduced number of 7 rays if the 
spectrum softens with energy. Because the imaging at- 
mospheric Cherenkov telescopes HESS and VERITAS 
have sensitivities of « 4 x 10~ 13 ph(> 300 GeV) cm~ 2 
s _1 for ~ 50 hour observations [ST] (see Fig. 2), M31, a 
northern hemisphere source (+41° declination), could be 
marginally detectable with VERITAS in long exposures, 
provided that the spectral GeV-TeV softening and re- 
duction in sensitivity due to M31's angular extent, ~ 1°, 
are not too great. 

The enhanced supernova rate in starburst galaxies 
such as M82 and NGC 253 at » 3 Mpc improve the 
prospects that they could be detectable with GLAST 
and ground-based Cherenkov telescopes |16j . In the inner 
starburst regions of these sources, interaction of cosmic 
rays with the strong stellar winds would produce GeV 
and TeV radiation [22]. These processes will also gener- 
ate neutrinos, though it is unlikely that they will be de- 
tected with IceCube or a Northern Hemisphere km-scale 
neutrino telescope, as these detectors have a sensitivity 
comparable to EGRET in terms of fluency- 
Torres (2004) [23] developed a detailed model of the 
nonthermal cosmic-ray production from the ultralumi- 
nous infrared galaxy (ULIRG) Arp 220 at « 72 Mpc. 
ULIRGs are the result of merging galaxies that drive 
large quantities of gas to the center of the system to 
form a dense gas disk, trigger a starburst, and possi- 
bly fuel a buried AGN. Because of their intense infrared 
emissions, Compton scattered radiations from cosmic ray 
electrons on the IR photons could additionally enhance 
the 7-ray fluxes. In spite of its large distance, Arp 220 
is potentially detectable with GLAST and the ground- 
based 7-ray telescopes [24], because the dense clouds of 
target gas and increased cosmic ray confinement signifi- 
cantly increase the brightness of ULIRGs in comparison 
to expectations from a simple scaling to the Milky Way. 

Clusters of Galaxies Nonthermal radiation will accom- 
pany the formation of collisionless shocks by merging 
clusters of galaxies during the merger of dark matter ha- 
los in the standard model for the growth of structure 
in a ylCDM universe [5S]. The available energy in the 
merger between a cluster of mass Mi = 10 15 Mi$Mq 
and a smaller cluster of mass M2 = 10 14 Mi4M Q , ini- 
tially separated by a distance of 7*1 = tm P c Mpc, is 



GM X M 2 8 x 10 1 



(;:>, 



r Mpc 



M 15 M 14 ergs 



(37) 




Fig. 2 Predicted 7-ray emission from the C oma cluster of 
galaxies from the cluster merger shock model |26| . The solid 
curve is the predicted photon flux (in units of ph(> -E pea k) 
cm -2 s _1 ). Sensitivity curves for EGRET, MAGIC, GLAST, 
and VERITAS and HESS are shown. The EGRET limits are 
for 2 weeks in the pointing mode, and the GLAST limits 
are for 1 year in its scanning mode. The quoted VERITAS, 
MAGIC and HESS point-source sensitivities are for 50 hour, 
5a observations 1211. 



If only a very tiny fraction of this energy is dissipated in 
the form of nonthermal cosmic rays protons, then £ 10 60 
ergs of cosmic rays, which are effectively trapped in the 
cluster even for a weak, ~ 0.1/xG, magnetic field, will 
be dissipated on the mean timescale for a nuclear colli- 
sion with cluster gas. The thermal X-ray bremsstrahlung 
cluster emission shows that the thermal cluster matter 
density is n t h sa 10~ 3 cm~ 3 , so that the characteristic 



timescale for nuclear interactions is t 
lis 



pp 



{n th a pp c) 



10 s. Hence the bolometric nonthermal cluster lumi- 



nosity from pion producing interactions is L pp 
-1 



10 4 



1 A fuller discussion of neutrino sources will be given in my 
Madison proceedings for TeV/Particle Astrophysics II. 



ergs s 

Fig. 2 shows Berrington's calculations [26 of the pre- 
dicted 7-ray emission from the Coma cluster of galax- 
ies, using parameters appropriate to the recent merger 
that has taken place in the Coma cluster environment 
(d = 100 Mpc). Point-source sensitivity limits for VER- 
ITAS (which is comparable to the HESS sensitivity) are 
taken from Ref. [21] . though the sensitivity may be de- 
graded by Coma's angular extent [28]. The predicted 7- 
ray emission falls below the EGRET sensitivity curve 
and the measured 2a upper limit of 3.81 x 10~ 8 ph(> 100 
MeV) cm~ 2 s^ 1 [27]. The results from the merging clus- 
ter model show that GLAST will significantly detect 
the non-thermal 7-rays from Coma to energies of several 
GeV. Furthermore, VERITAS could have a high confi- 
dence ( £ 5fi) detection of Coma (at declination +38°), 
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Fig. 3 Diffuse extragalactic 7-ray background, from analyses 
of EGRET data [19U20] . compared to model calculations of 
the contributions to the EGRB for FSRQ and BL blazars, and 
tot al A GNs [5], star- forming galaxies [32J . star burs t galax- 
ies [33], structure shocks in clusters of galaxies [30] . and all 
long-duration GRBs, including those detected as such (esti- 
mated herein). Pulsar contribution at 1 GeV is ~ 20% of 
star-forming galaxy estimates. 

depending on detail on the fraction of energy going into 
shocked cosmic-ray protons, the nonthermal nature of 
Coma's hard X-ray spectrum, and the amount of non- 
thermal proton energy left over from previous merger 
events. 

In addition to merger shocks, high Mach number ac- 
cretion shocks at the periphery of a forming cluster can 
accelerate nonthermal particles [29]. The diffuse 7-ray 
background formed by intergalactic structure formation 
shocks from the calculations of Keshet et al. (2003) [30] 
is shown in Fig. 3. Gamma-ray emission has not yet been 
convincingly detected from clusters of galaxies |27j . but 
calculations like this indicate that clusters of galaxies 
are likely to be the next established class of extragalac- 
tic sources of high-energy radiation. 



Using the star formation rate function of Ref. [7] , the 
integral in eq. ([3"5| is easily performed to give a value of 
2.14, so that el* 8.7 x 10~ 7 [£: 7 /(100 MeV)]-° 4 GeV 
cm" 2 s" 1 sr" 1 . This result is in good agreement with the 
more detailed treatment of the "guaranteed 7-ray back- 
ground" by Pavlidou and Fields (2002) [35], plotted in 
Fig. 3. The upper curve is scaled to a dust-corrected star 
formation rate, and the lower curve stops the integration 
at z = 1 . Also shown is the intensity of starburst galaxies 
estimated from a radio/FIR correlation [3"5] . 

Extragalactic Pulsar Emissions Using SAS-2 data for 
the diffuse galactic 7-ray emission and EGRET data for 
pulsars, the analysis of Ref. pH] shows that the total 
> 100 MeV flux of diffuse radiation from the Milky Way 
is w 1.5 x 10~ 7 ergs cm -1 s" 1 , and the combined flux 
of the 6 brightest EGRET pulsars is » 1.35 X 10~ 8 ergs 
cm -1 s" 1 . The modeling in that paper shows that the 
superposition of diffuse fluxes of unresolved pulsars is at 
the level of ss 1.2 x 10~ 8 ergs cm -1 s" 1 . Thus total pulsar 
emissions make up as much as 20% of the total galactic 
7-ray flux in star-forming galaxies like the Milky Way. 

The apparently diffuse emissions from pulsars in galax- 
ies throughout the universe is then, to first order, at the 
level of w 20% and proportional to the star formation 
history of the universe. A number of important effects 
must be considered for more accurate estimates of the ex- 
tragalactic diffuse pulsar flux at different 7-ray energies, 
most obviously being the harder pulsar spectrum (com- 
pared to the cosmic-ray induced emissions) at energies 
up to the pulsar cutoff energies between « 1 - 100 GeV 
[3"51l36| . Of great interest is to accurately measure the 
high-energy pulsar spectral cutoffs with GLAST, which 
can be included in a more complete model for the pulsar 
contribution to the diffuse galactic background. Unfortu- 
nately, GLAST would not be sensitive to detect Milky- 
Way like 7-ray pulsars from nearby galaxies. Placing the 
Crab pulsar at 1 Mpc would yield a ^ 100 MeV appar- 
ent isotropic flux ,$ 10~ 14 ergs cm" 2 s" 1 . 



Diffuse Intensity from Star-Forming Galaxies We can 
use eq. (|28p to calculate the diffuse intensity from star- 
forming galaxies by normalizing to the density and non- 
thermal 7-ray luminosity of L* galaxies like the Milky 
Way. Fits of galaxy surveys to the Schechter luminos- 
ity function imply that the density of L* galaxies is 
n* = 0.016ft 3 Mpc" 3 « 1/(170 Mpc 3 ) [31]. Employing a 
mono-luminosity galaxy luminosity function, e^q^e*; z) 
= n*£*(z)e*L*(e*), and eq. ([28]) becomes, using eq. ([34]) , 



ei: = — 



c 

-in 



dz — 

dz 



lq*(e*;z) 
l + z 



4x10"' x 



V 200^ J 



dz 



S,{z)(l + z)- 2A 

y/n m (l + Zf + f2 A 



GeV 



cm' s sr 



(38) 



5.2 7-rays and v from Beamed Sources 

The evidence from EGRET and Whipple shows that 
beamed GRBs and blazars are the brightest extragalactic 
high-energy 7-ray sources, and that isotropically emit- 
ting sources will be difficult to detect except in a few 
cases, as just demonstrated. 



Microquasars One class of beamed source that has not 
yet been detected from beyond the Galaxy is the mi- 
croquasar class, even though some of the ultraluminous 
X-ray sources seen in nearby galaxies could be micro- 
quasars with their jets oriented towards us [37] . Presently 
only high-mass microquasars are known sources of GeV 
and TeV radiation, and the evidence of associations of 
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7-ray sources with low-mass systems is weak. The estab- 
lished cases of 7-ray emitting microquasars are LS 5039, 
associated with an unidentified EGRET source [38] and 
unambiguously detected with HESS [35], and LSI +61 
303, whose orbital modulation has been recently demon- 
strated with the MAGIC telescope [40] . 

Models for microquasars do not require large bulk 
Lorentz factors r of the plasma outflow in microquasar 
jets [41j . and superluminal motion observations gener- 
ally reveal microquasar Doppler factors 2. The small 
sample leaves open the possibility that r could exceed 
a few, which would make possible the detection of mi- 
croquasars from distant galaxies. Using eq. (TIT)]) and the 
threshold expression 8^' e t aiJ / d 2 f e = 10~ 12 ergs cm~ 2 
s _1 , we find that for a flat vF v spectrum (a v = 2) that 
tip }t 6(d/Mpc) 2 for a Chandrasekhar mass compact ob- 
ject and a bolometric factor Xf, = 10. Thus, Doppler fac- 
tors of only a few or greater are needed in order to detect 
microquasars from galaxies at distances ^ 1 Mpc. 

Gamma Ray Bursts A detailed treatment of the statis- 
tics of GRBs is given in Ref. [5J. We can use those re- 
sults and eq. ([25]) to estimate the diffuse 7-ray back- 
ground from GRBs. To simplify the results, we consider 
all GRBs, including those above threshold, and approx- 
imate the GRB spectrum as a flat vF v spectrum to the 
highest 7-ray energies. The diffuse 7-ray background in- 
tensity is then given by 

47r\ b H J dZ (1 + z) 2 ■ m 
For a local GRB event rate tigrb — lOriio Gpc~ 3 yr , 
grb ~ 3 x 10~Vo 



el 



GRB r^j al GRB^* 



el 



dz aai, (40) 



A10 Ja ^jQ m {\ + zf + Q A 

in units of GeV cm -2 s _1 sr -1 . 

Using star formation rates 5 and 6 that allow the 
redshift and jet opening angle distributions to be fit [BJ 
gives the diffuse intensity from GRBs shown in Fig. 3 
which is a small fraction of the diffuse 7-ray background. 

Blazars EGRET and GLAST data on blazars can be an- 
alyzed with models that jointly fit the redshift and flux 
size distribution and predict the level of the extragalac- 
tic 7-ray background (EGRB) 0H]. The diffuse intensi- 
ties of unresolved FSRQ (dashed) and BL Lac (dotted) 
blazars, and the total AGN contribution (shaded), are 
shown in Fig. 3 [5]. The ranges correspond to sensitivi- 
ties <f>-$ — 25 and 0_g = 12.5. The BL Lac objects and 
FSRQs, including emissions from misaligned radio galax- 
ies, contribute at the ~ 2 - 4% and ~ 10 - 15% levels, 
respectively, to the total EGRB [1]5] near 1 GeV. 

The sum of the different contributions in Fig. 3 at 
~ 1 GeV is at about the level of the total diffuse ex- 
tragalactic 7-ray emissions measured with EGRET [19] . 
Soft blazar sources, and softer than modeled diffuse cos- 
mic ray emissions from normal galaxies, could account 



for the residual emissions between w 50 MeV - 1 GeV. 
New hard 7-ray source populations are apparently re- 
quired at £ 10 GeV, which would include cascade emis- 
sion from UHE electromagnetic cascades with photons 
of the EBL. 



6 Summary and Conclusions 

It will be of considerable interest when GLAST or a 
TeV telescope detects M31 or another galaxy of the local 
group, or a starburst or IR luminous galaxy. The mea- 
sured flux will give a valuable check on the efficiency of 
cosmic ray production as a function of galactic star for- 
mation activity, and will provide a further normalization, 
after the Milky Way and the LMC, of the contribution 
of star forming galaxies to the 7-ray background. 

GLAST will provide large statistical samples on at 
least two source classes: GRBs and blazars. There are 
good reasons to think that GLAST will detect star form- 
ing galaxies and clusters of galaxies. Although the sensi- 
tivity to extended sources is degraded in TeV telescopes, 
the Coma cluster is near the threshold for detection, 
depending on its nonthermal X-ray spectrum, but M31 
would require fortunate spectral and spatial emissions to 
be detected. 

For GRBs and blazars, there are already fundamen- 
tally interesting questions about the evolution of the 
source rate densities and change in properties of relativis- 
tic jet sources through cosmic time. For 7-ray emitting 
BL Lac objects, positive source evolution (more sources 
at late times) and negative luminosity evolution (sources 
brighter in the past) explains the statistical distributions 
from EGRET [5]. Detection of blazars to high redshift, 
z ^> 5, is expected with GLAST. It will be interesting to 
see if GLAST LAT-detected GRBs are peculiar in their 
properties compared to long-duration GRBs found by 
burst detectors with ~ 100 keV triggers. 

The EGRB is probably a composite of many source 
classes (Fig. 3), which can best be established by identi- 
fying individual sources with better sensitivity detectors. 
The contribution of beamed to unbeamed sources, after 
subtracting known sources, is limited by the 7-ray sta- 
tistical excursions of the EGRB measured with GLAST. 
GLAST will monitor blazar and GRB flaring with its 
~ 2 sr field of view, and rapidly slewing instruments like 
MAGIC may soon discover the first VHE, > 10 GeV, 
GRB. But a wide field-of-view ground-based 7-ray tele- 
scope, like HAWC or the next generation TeV telescopes, 
will have the best chance to monitor TeV 7-ray transients 
for follow-up observations. Improved statistical analyses 
of 7-ray and particle astronomy projects will tell us the 
composition of the unresolved residual 7-ray emission, 
and if there are room for more source classes, such as 
dark matter emissions, anomalous microquasars or odd 
classes of 7-ray emitting objects yet to be discovered. 
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